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2 Abstract 

We consider the questions of whether the damped Lyman-alpha (DLA) and 
sub-DLA absorbers in quasar spectra differ intrinsically in metallicity, and 

t-h whether they could arise in galaxies of different masses. Using the recent 

measurements of the robust metallicity indicators Zn and S in DLAs and 
sub-DLAs, we confirm that sub-DLAs have higher mean metallicities than 
DLAs, especially at z < 2. We find that the intercept of the metallicity- 
redshift relation derived from Zn and S is higher than that derived from 

q Fe by 0.5-0.6 dex. We also show that, while there is a correlation between 

the metallicity and the rest equivalent width of Mg II A2796 or Fe II A2599 
for DLAs, no correlation is seen for sub-DLAs. Given this, and the similar 
Mg II or Fe II selection criteria employed in the discovery of both types 
of systems at lower redshifts, the difference between metallicities of DLAs 
and sub-DLAs appears to be real and not an artefact of selection. This 
conclusion is supported by our simulations of Mg II A2796 and Fe II A2599 
lines for a wide range of physical conditions. On examining the velocity 
spreads of the absorbers, we find that sub-DLAs show somewhat higher mean 
and median velocity spreads (A-y), and an excess of systems with Av > 150 
km s -1 , than DLAs. Compared to DLAs, the [Mn/Fe] vs. [Zn/H] trend for 
sub-DLAs appears to be steeper and closer to the trend for Galactic bulge 
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and thick disk stars, possibly suggesting different stellar populations. The 
absorber data appear to be consistent with galaxy down-sizing. The data are 
also consistent with the relative number densities of low-mass and high-mass 
galaxies. It is thus plausible that sub-DLAs arise in more massive galaxies 
on average than DLAs. 

Keywords: (galaxies:) quasars: absorption lines-galaxies: evolution-ISM: 
abundances 



1. Introduction 

The evolutionary history of galaxies is written in part in the interstellar 
abundances of the chemical elements. Deciphering this history from observa- 
tions of distant galaxies is an important aspect of understanding the various 
processes responsible for evolution of galaxies. Besides the stellar mass func- 
tion, yields, supernova rates etc., the metallicity of a galaxy is also governed 
by the depth of its gravitational potential well, which determines the ability 
to retain the metals produced from star formation. Indeed, the interstel- 
lar metallicity (as determined from nebular emission lines) is found to be 
correlated tightly with the stellar mass of the galaxy at z ~ and also at 
0.4 < z < 3 (e.g., Tremonti et al. 2004; Savaglio et al. 2005; Erb et al. 2006). 

Quasar absorption lines provide a powerful complementary tool for study- 
ing distant galaxies. The damped Lyman-alpha (DLA; log N m > 20.3) and 
sub-damped Lyman-alpha (sub-DLA; 19.0 < log N m < 20.3) absorbers are 
the primary neutral gas reservoir at < z < 5 (e.g., Storrie-Lombardi & 
Wolfe 2000; Peroux et al. 2005; Prochaska et al. 2005), and offer the most 
precise element abundance measurements in distant galaxies. 

As per most chemical evolution models, the mean interstellar metallicity 
of galaxies should reach a near-solar value at low redshift, as a result of pro- 
gressive generations of star formation. Surprisingly, most DLAs observed at 
0.1 < z < 3.9 are metal-poor (e.g., Kulkarni et al 2005, 2007, and references 
therein). By contrast, a substantial fraction of the sub-DLAs observed at 
0.6 < z < 1.5 are metal-rich-some even super-solar (e.g., Kulkarni et al. 
2007; Meiring et al. 2007, 2008, 2009a, 2009b; Peroux et al. 2006, 2008; 
and references therein). Based on the observed mass-metallicity relation for 
galaxies, Khare et al. (2007) suggested that the sub-DLAs may be more 
massive compared to DLAs. Here we examine various selection effects that 
could potentially be responsible for this conclusion with an updated sample 
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of DLA/sub-DLA metallicities, and compare other properties of DLAs and 
sub-DLAs. We use the nearly undepleted elements Zn and S, which makes 
our analysis "cleaner" and more conservative than studies based on depleted 
elements such as Fe. This also allows us to make meaningful comparisons 
of the observed DLA/sub-DLA metallicity evolution with known trends for 
galaxies. 

2. Sample Description, Observed Trends, and Metallicity Indicator 

2.1. Sample Description 

Comparison of the metallicity evolution of DLA and sub-DLA galaxies 
was made in Kulkarni et al. (2007) based on Zn data for 119 DLAs and 
30 sub-DLAs. In this paper, we examine this issue further using Zn and S 
for 154 DLAs and 58 sub-DLAs. The two-fold increase in the number of 
sub-DLAs since 2007 has come primarily from our recent measurements (20 
systems from Meiring et al. 2008; Peroux et al. 2008; Meiring et al. 2009a), 
complemented by a few measurements from the literature (Nestor et al. 2008, 
Quast et al. 2008, Noterdaeme et al. 2008, Dessauges-Zavadsky et al. 2009). 
To avoid any prior bias toward high or low metallicity, we exclude systems 
which were specifically observed because of hints that they are rich in metals 
or molecules. (Thus, we exclude two metal-rich sub-DLAs from Prochaska et 
al. 2006, one metal-rich sub-DLA from Srianand et al. 2008, and the metal- 
strong DLAs from Herbert-Fort et al. 2006. In any case, we note Kulkarni 
et al. 2007 showed that adding the metal-strong systems did not change the 
conclusions about relative evolutions of DLAs and sub-DLAs.) 

2.2. Observed Trends 

Fig. 1 shows the iVm-weighted mean metallicity vs. look-back time rela- 
tion for the 154 DLAs and 58 sub-DLAs in our sample, calculated using the 
procedures outlined in Kulkarni & Fall (2002). Upper limits on Zn have been 
treated with survival analysis. Also shown for reference are the predictions for 
metallicity evolution in two theoretical models. The light dot-dashed curve 
shows the mean interstellar metallicity in the chemical evolution model of 
Pei et al. (1999) with the optimum fit for the cosmic infrared background 
intensity. This model also uses observational constraints from optical galaxy 
surveys and the comoving H I density of DLAs to calculate the coupled 
global evolution of stellar, gaseous, and metal content of galaxies. The light 
dot-double-dashed curve shows the mean metallicity of cold interstellar gas 
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in a semianalytic model of galaxy formation in the cold dark matter merg- 
ing hierarchy by Somerville et al. (2001). This model, referred to as the 
"collisional starburst model" by Somerville et al., assumes star formation in 
bursts triggered by galaxy mergers in addition to a quiescent star formation 
at constant efficiency 

The sub-DLA global mean metallicity appears to be higher than that of 
DLAs, reaching a near-solar value at low z, consistent with the models (see 
also York et al. 2006; Prochaska et al. 2006). The linear regression slope 
for the sub-DLA Ani-weighted mean metallicity vs. redshift data (—0.46 ± 
0.18) appears to be steeper than that for DLAs (—0.22 ± 0.07). The slope 
for sub-DLAs has larger uncertainties due to the smaller sample size. The 
corresponding linear regression estimates of the intercepts, i.e. expected 
metallicities at z — 0, are 0.18 ± 0.29 for sub-DLAs and —0.63 ± 0.16 for 
DLAs, which differ at 2.4 a level. The bold curves in Fig. 1 show the best 
fits to the DLA and sub-DLA data (plotted in terms of look-back time). 
The results do not change substantially if the metallicities are unweighted 
(a slope of —0.47 ± 0.16 and an intercept of 0.22 ± 0.24 for sub-DLAs, and 
a slope of —0.25 ± 0.05, intercept of —0.47 ± 0.11 for DLAs.) Below some 
redshift around 2, the mean metallicity of sub-DLAs exceeds that of DLAs. 
Determining exactly where that happens (e.g. at z < 1.7 or z < 2.2) will 
require a larger database. 

2. 3. The Choice of the Metallicity Indicator 

We have used measurements of Zn or S, since these nearly undepleted 
elements offer the most direct "dust-free" metallicity indicators. Further- 
more S and Zn are seen to track each other well (within < 0.2 dex) in most 
absorbers, as seen in Fig. 2, which shows [S/Zn] vs. [Zn/H] for absorbers 
with detections of both S and Zn (see also Nissen et al. 2004, 2007). Fe, 
used in some previous works, is not a reliable metallicity indicator since it is 
severely depleted on interstellar dust, so that its gas-phase abundance does 
not give its total abundance (see, e.g., Miller et al. 2007; Jenkins 2009). 
Moreover, depletions correlate with metallicity: [Zn/Fe] ~ 0.2 for [Zn/H] 
~ — 1.5 and [Zn/Fe] ~ 1 for [Zn/H] ~ and can show scatter (e.g., Ledoux 
et al. 2003; Meiring et al. 2006). Thus, [Fe/H] can seriously underestimate 
the true metallicity, especially at the high-metallicity end, and may suggest 
a less steep evolution than the true trend. It is interesting that the slope for 
sub-DLAs from unweighted Fe/H (—0.41 ± 0.07; Dessagues-Zavadsky et al. 
2009) is close to our value —0.47 ± 0.16 based on unweighted X/H (where 
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X = Zn or S), but the intercepts based on unweighted Fe/H (—0.26 ± 0.13 
and —1.03 ± 0.09 for sub-DLAs and DLAs, respectively) are substantially 
lower than those based on unweighted X/H (0.22 ±0.24 and — 0.47±0.11 for 
sub-DLAs and DLAs, respectively). Thus, using Fe systematically underes- 
timates the intercept of the metallicity-redshift trend by 0.5-0.6 dex for both 
DLAs and sub-DLAs. 

3. Is the Metallicity Difference between DLAs and Sub-DLAs a 
Selection Effect? 

We now consider whether the higher metallicity of sub-DLAs could be 
an artefact. Ionisation effects and dust selection have been ruled out before 
as the cause of the difference between DLA and sub-DLA metallicities (e.g., 
Khare et al. 2007). We mainly consider here the possiblity of a selection 
effect tied to the presence of strong Mg II A 2796 and Fe II A2599 lines. 

3.1. Is there an Mg II Selection Bias? 

The sample of sub-DLAs studied by our team over the past few years 
(Meiring et al. 2009b and references therein) was chosen from the low-z 
sub-DLAs with measured Nm values. As emphasized in our previous papers, 
the availability of JVhi was the primary selection criterion, not the strength 
of any metal lines. However, the HST sample of DLAs and sub-DLAs that 
S. Rao and collaborators chose to measure Nm are systems with rest-frame 
equivalent widths for the Mg II A2796 line exceeding 0.5 A (14^796 > 0.5 
A). Below we discuss whether this introduces a bias toward more metal-rich 
systems in the sub-DLAs. Simply put, the argument is that a system with 
large W2796 would have large A^Mgii, hence a large iVz n n, and hence a higher 
metallicity for low Nm- There are several problems with this argument: 

Mg II Line Saturation: The Mg II A2796 line is almost always highly 
saturated in DLAs and sub-DLAs. The systems with larger W^2796 have larger 
velocity spreads, not larger Mg II column densities. Thus, selecting systems 
with larger W2796 does not necessarily translate into selecting systems with 
larger N MgU and hence larger N Z nu- 

Lack of Correlation between Mg II Strength and Metallicity: 
If there were an Mg II selection bias, one would expect to see a positive 
correlation between metallicity and H2796 for sub-DLAs. However, no corre- 
lation is seen between [X/H] (where X = Zn or S) and W 2 7ge f° r sub-DLAs 
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(see Fig. 3a). Using the generalized Kendall's r-correlation coefficient (in- 
cluding survival analysis to treat upper limits), we get r = —0.070, and a 
2-sided probability of 0.670 for no correlation (see Table 1). Even [Fe/H] and 
W2796 do not correlate for sub-DLAs, with generalized Kendall's r = 0.128 
and a 2-sided probability of 0.541 for no correlation for points from Fig. 3 
of Dessagues-Zavadsky et al. (2009). By contrast, a significant correlation 
between Mg II strength and [X/H] is seen for DLAs (Fig. 3b), with a gener- 
alized Kendall's r = 0.819, and a probability of < 0.0001 for no correlation 
(Table 1). Even [Fe/H] shows a significant correlation with \V279e for DLAs 
(see also Fig. 1 of Murphy et al. 2007). Thus, the sub-DLAs appear to be 
far less likely than the DLAs to show an Mg Il-selection bias toward higher 
metallicity. The smaller correlation in the sub-DLA sample is apparent even 
looking at just the detections plotted in Figs. 3(a) and 3(b). It is possible 
that this difference between DLAs and sub-DLAs is caused partly by the 
smaller number of sub-DLAs plotted in Fig. 3(a), but the number of DLAs 
in Fig. 3(b) is not much larger, and the DLAs appear to show a correlation 
across most of the range of PU2796- 

Lack of Sensitivity Problem: Kulkarni et al. (2007) used simulations 
of the Mg II A2796 line with a variety of velocity component structures to 
estimate the minimum measurable [Mg/H] in the observed sub-DLA sample. 
These simulations showed that, for the choice of W2796 > 0.5 A one can 
detect systems with [Mg/H] well below —1.5 for log N m = 19.7. Thus, 
choosing systems with W279Q > 0.5 A cannot explain the fact that ~ 50% of 
the \ow-z sub-DLAs with Zn or S detections have [X/H] > —0.5, and only 
~ 10% systems with Zn detections have [X/H] < —1.0. 

Representativeness of Sub-DLA Sample: One may wonder whether 
the sub-DLAs with > 0.5 A are unrepresentative of the sub-DLA 

population in terms of metallicity. However, for this to produce a bias in the 
current measurements, the number distribution of sub-DLAs would have to 
be skewed to have many more systems with W2796 < 0.5 A than with W2796 > 
0.5 A, and those objects would have to be more metal-poor than what is 
found now for sub-DLAs. This does not seem probable. More important, 
such a scenario would imply a large missing fraction of sub-DLAs that would 
be difficult to reconcile with current measurements of the number density 
dn/dz and the column density distribution f(Nm) which are consistent with 
both DLAs on the one hand and Lyman limit systems on the other hand. 

Metallicity correlation with E(B — V): York et al. (2006, hereafter 
Y06), using > 800 SDSS Mg II systems and assuming a constant dust-to- 
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gas ratio, found an anti-correlation between quasar reddening E{B — V) and 
metallicity. This trend becomes more pronounced for a metallicity-dependent 
dust-to-gas ratio (Khare et al. 2007). Y06 also showed that strong \V279e can 
occur for all values of E(B — V), which implies that systems with strong 
Mg II lines could correspond to high or low metallicity, i.e. using systems 
with large W 2 im does not bias one toward selecting only metal-rich systems. 

We thus conclude that the higher sub-DLA metallicities are not caused 
by an Mg II selection bias. 

3.2. Is there an Fe II Selection Bias? 

Another potential bias could come from the fact that the \ow-z DLA 
and sub-DLA absorbers chosen for the HST UV spectroscopy by Rao et al. 
also have relatively strong Fe II A2599 lines (rest frame equivalent widths 
W 2599 > 0.5 A). The squares and circles in Figs. 4(a) and 4(b) show the 
observed metallicity vs. 1^2599 for the sub-DLAs and DLAs, respectively, 
in our sample for which IF2599 measurements are available. A correlation 
between the observed metallicity and the strength of Fe II A2599 is seen for 
DLAs, but not for sub-DLAs (Table 2). This shows that the sub-DLA data 
have no systematic bias toward high metallicity systems. 

To examine this issue further, we explored the variation of ^2599 with 
metallicity by simulating Fe II A2599 lines in DLAs and sub-DLAs consid- 
ering a range of physical conditions of the absorbing gas. The H I column 
densities were chosen to resemble the ranges observed in our sample: for 
DLAs, randomly in the range 20.3-22.0, and for sub-DLAs, such that 3/4 of 
all the randomly selected values lie in the range 19.7-20.3, and 1/4 in the 
range 19.0-19.7. For each selected value of N ffl , the metallicity [X/H] was 
selected randomly between -2.5 and +1.0, and the depletion factor for Fe 
was selected randomly within a range of -0.6 to -1.2 dex. [This is a fairly 
conservative Fe depletion range in warm Galactic ISM. The \ow-z sub-DLAs 
with Zn detections show [Fe/Zn] in the range of -0.1 to -1.6 dex, most be- 
tween -0.4 and -0.9 dex (e.g., Meiring et al. 2009b).] The number of velocity 
components in each artificial profile was chosen randomly between 3 and 15 
having random velocity values within a range constrained by the selected 
number of components (the maximum range being -375 to 375 km s _1 for 15 
components). The components were assigned random Doppler b parameters 
ranging between 3 and 20 km s _1 . The equivalent width of each resulting 
Fe II A2599 line was measured. Finally, the results for the 30,000 artificial 
lines generated in this way were binned over metallicity into 0.25 dex bins. 
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The diamonds in Figs. 4(a) and 4(b) show the median equivalent width in 
each bin plotted vs. the central metallicity of the bin. The error bars on the 
median equivalent widths in each metallicity bin show the ±ler variations 
among the equivalent widths generated from the simulations with metallici- 
ties lying in that bin. Over ~ 2/3 of the observed detections lie within the 
±lcx error bars, as may be expected. 

Figs. 4(a) and 4(b) suggest that one may expect a gentle trend of gradu- 
ally increasing W2599 with increasing metallicity. Indeed, the observed DLAs 
do show a correlation roughly resembling the median simulated trend, except 
that no supersolar DLAs are seen. By contrast, although most sub-DLA 
detections lie within the ±la range of the simulated trend, they show no re- 
semblance to the median simulated trend. The lack of correlation among the 
sub-DLA points may be a result of the small sub-DLA sample. The cluster 
of points lying to the right of the simulated trend in Fig. 4(a) suggest that 
the depletions in some of the observed sub-DLAs may be less severe than 
assumed above, or the velocity components may be more widely spread. In 
any case, it is clear from Fig. 4(a) that a sample with IF2599 > 0.5 A can 
have sub-DLAs with [X/H] as low as -1.1 even for the relatively large deple- 
tion factors assumed above (as low as -1.3 within the 2 a range). All of the 
observed metallicities for sub-DLAs with detections are considerably higher 
than that. Indeed, Fe II selection does not explain why several sub-DLAs 
have super-solar metallicity, while DLAs do not. 

We conclude that the higher mean sub-DLA metallicities compared to 
DLAs do not appear to arise from an Fe II selection bias. This conclusion 
is consistent with the lack of Mg II selection bias (section 3.1), since Fe II A 
2599 is fairly close in strength to Mg II A 2796 (log [A x /] = 2.793 for Fe II 
A 2599 and 3.236 for Mg II A 2796). 

3.3. Do abundance gradients play a role? 

It has been suggested that DLAs could be metal-poor because gas cross- 
sections are larger in the less developed outskirts of galaxies (e.g., Chen et al. 
2005; Ellison et al. 2005; Zwaan et al. 2005). However, abundance gradients 
are relatively weak in the Galactic disk, e.g., an [Fe/H] gradient of -0.052 
dex kpc -1 for 5-17 kpc and -0.012 dex kpc -1 for 10-15 kpc (Lemasle et al. 
2008). The SDSS has shown that the stellar metalicity gradient at the solar 
circle is low, < 0.01 dex kpc -1 (Ivezic et al. 2008). Metallicity gradients are 
also weak in nearby spirals, ~ -0.09 dex kpc -1 (Ferguson et al. 1998), and 
it is unclear whether the gradients exist in distant galaxies. In fact, given 
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that H I column densities are lower in sight-lines through outer parts of local 
galaxies, one may expect sub-DLAs to have lower abundances than DLAs. 
Thus abundance gradients cannot explain the large difference between the 
metallicities of DLAs and sub-DLAs. 

How well-determined are the H I column densities? As discussed by 
Meiring et al. (2007, 2008, 2009a-see appendices) and by Ellison, Murphy, 
& Dessauges-Zavadsky et al. (2009a), the Nm values of sub-DLAs inferred 
from single-component fits to low-resolution spectra can be overestimates 
by 0.1-0.3 dex. (If there is a wide spread of components and Ahi is low, a 
single-component H I fit forces the damping profile to contain more H I that 
is actually there. The effect is not important for DLAs.) Thus, the true 
metallicities of sub-DLAs can be higher by 0.1-0.3 dex, making the difference 
between DLAs and sub-DLAs even larger. 

We conclude that sub-DLAs are, on average, intrinsically more metal-rich 
than DLAs. Near-solar or super-solar abundances have also been estimated 
in some Lyman-limit systems and even some systems with log Am < 16 
(e.g., Ding et al. 2003; Maseiro et al. 2005; Jenkins et al. 2005; Misawa 
et al. 2008). These studies suggest (despite the dependence on details of 
photoionisation calculations) the intriguing possibility of the Ani-metallicity 
anti-correlation extending over ~ 6 orders of magnitude in Ahi- 

4. The Velocity Spreads of DLAs and Sub-DLAs 

One may wonder whether the metallicity difference between DLAs and 
sub-DLAs is related to kinematic differences. We therefore compare the 
velocity spreads of the absorbers. 

4-1- Velocity Spread Distributions 

The velocity spread Av for DLAs and sub-DLAs is usually derived by 
integrating the apparent optical depth profile over the extent of a metal line, 
and calculating the velocity width enclosing the inner 90% of the absorption. 
While the overall range of Av is similar for our DLA and sub-DLA samples, 
the mean velocity spread is 125.5 km s" 1 for sub-DLAs and 103.0 km s _1 
for DLAs. The median values are 95.3 km s _1 for sub-DLAs and 80.0 km 
s _1 for DLAs. Thus sub-DLAs have somewhat larger velocity spreads on 
average than DLAs. Furthermore, the fraction of systems with Av > 150 
km s" 1 relative to those with Av < 50 km s" 1 is greater among sub-DLAs 
than DLAs (see Table 3). Assuming Poissonian error bars on the number 
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of systems detected in either velocity range, sub-DLAs show a 2.0 a excess 
of systems with Av > 150 km s _1 as compared to DLAs. Of course, larger 
samples of both DLAs and sub-DLAs are needed to verify the difference. 

4-2. Velocity-Metallicity Correlations 

Correlations between velocity spreads and metallicities might exist for 
physical reasons. For instance, lower mass galaxies might develop higher wind 
velocities in outflows. Ledoux et al. (2006) reported a velocity-metallicity 
correlation in a sample of 62 DLAs combined with 8 high-iVm sub-DLAs, 
and suggested that it reflects a mass-metallicity relation. Figs. 5(a) and 
5(b) show velocity spread vs. metallicity for sub-DLAs and DLAs in our 
sample that have velocity measurements available. A generalized Kendall's 
t correlation test reveals that the correlation of [X/H] with Av is significant 
for DLAs, but not for sub-DLAs (Table 4). We find that the correlation 
between [Fe/H] and Av from Fig. 2 of Dessagues-Zavadsky et al. (2009) is a 
little stronger than the trend between [X/H] and Av, but still marginal for 
sub-DLAs. A weaker correlation with Av for [X/H] than for [Fe/H] could be 
partly due to the lack of a significant number of systems with low [Zn/H], 
or it could be caused by the Routly-Spitzer effect (Routly & Spitzer 1952): 
systems with larger Av may have more supernova-driven shocks capable of 
destroying the interstellar dust grains, resulting in lower depletion of iron, 
and hence higher [Fe/H]. 

We note that there is no unique or perfect way to measure the velocity 
width, which is very sensitive to the S/N and the spectral resolution. Fur- 
thermore, for both DLAs and sub-DLAs, the mean or median Av derived 
from Mg II A2796 are about a factor of 2 larger than those derived from the 
weaker lines (Table 3). As an extreme example, we note that the sub-DLA 
toward Q05 15-4414 has Av = 573 km s" 1 from the widely spread Mg II A 
2796 line, but only 94 km s^ 1 from the weaker lines! Indeed, if Av from Mg 
II A2796 is used, the correlation with metallicity [X/H] becomes less robust 
for DLAs as well as sub-DLAs (see Table 4), although this could be partly 
due to the much smaller number of systems with [X/H] that also have Av 
measured from Mg II lines. Thus we caution against strong claims and in- 
terpretations related to the velocity spreads of DLAs and sub-DLAs. But we 
note that the larger velocity spread in Mg II A2796 indicates the presence of 
high-velocity components that are hard to see in the weaker lines, possibly 
suggesting a substantial amount of halo gas or high velocity clouds. (Mg II 
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in the high velocity clouds in the Milky Way can be inferred from other low 
ions such as Ca II and Si II [e.g., Wakker et al. 2007; Shull et al. 2009].) 

5. Can Metal-rich Sub-DLAs be more Massive than Metal-poor 
DLAs?: Comparison with Galaxy Data 

Based on the higher mean metallicity of sub-DLAs and the mass-metallicity 
relation for galaxies, Khare et al. (2007) suggested that sub-DLAs may arise 
in more massive galaxies, on average, than DLAs. This suggestion is not that 
shocking, once one realizes that a higher H I column density by no means im- 
plies a higher total mass. Indeed, the suggestion seems reasonable given that 
sub-DLAs have somewhat higher velocity spreads on average than DLAs, 
and that a higher fraction of sub-DLAs have Av > 150 km s _1 (see Table 
3). (Of course, as stressed by Kulkarni et al. 2007, not every sub-DLA is 
metal-rich. Thus, it is not necessary that every sub-DLA should be massive 
and have a large velocity width.) 

To further explore whether sub-DLAs can arise in more massive galaxies 
than DLAs, we now compare the metallicity evolution of DLAs and sub- 
DLAs with galaxy observations. Our use of the robust metallicity indicators 
Zn and S enables us to make meaningful comparisons with galaxy data. We 
show that the possibility of sub-DLAs being more massive than DLAs is 
consistent with stellar [Mn/Fe] ratios in galaxies, with galaxy down-sizing, 
with evolution of luminous star forming galaxies, and with galaxy number 
densities. 

Comparison with [Mn/Fe] ratios in galaxies: In the stars of the 
Milky Way, [Mn/Fe] is found to increase with increasing [Fe/H], which could 
indicate production of Mn in type la supernovae (e.g., Samland 1998), or 
metallicity-dependent yields in type II supernovae (e.g., Timmes, Woosley, 
& Weaver 1995; McWilliam et al 2003; Feltzing et al. 2007). Mc William et 
al. (2003) found that the [Mn/Fe] vs. [Fe/H] trend for the Sagittarius dwarf 
spheroidal is lower than that for the Galactic bulge stars by ~ 0.2 dex. They 
interpreted this difference to be a result of less rapid chemical evolution of 
the Sagittarius dwarf compared to the Galactic bulge: the Galactic bulge 
experienced rapid enrichment with time scales of < 1 Gyr (e.g., Matteucci, 
Romano, & Molaro 1999), while dwarf spheroidal galaxies have extended star 
formation histories lasting many Gyrs (e.g, Mateo 1998; Grebel 2000). 

DLAs and sub-DLAs also show increasing [Mn/Fe] with increasing [Zn/H]. 
However, the two trends differ, suggesting different stellar populations. Based 
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on our measurements from Meiring et al. (2009a), we find that the lin- 
ear regression slopes of the [Mn/Fe] vs. [Zn/H] trends are 0.27 ± 0.03 and 
0.06 ± 0.03, respectively, for sub-DLAs and DLAs. The slopes of the corre- 
sponding relations for the Galactic bulge stars in McWilliam et al. (2003) and 
the thick disk stars in Reddy et al. (2006) are 0.40 ± 0.05 and 0.34 ± 0.02 re- 
spectively. The slope for the stars in the Sagittarius dwarf spheroidal galaxy 
in McWilliam et al. (2003) is 0.15 ± 0.04. The [Mn/Fe] vs. [Zn/H] trend 
in sub-DLAs appears to be closer to that for bulge stars of the Milky Way, 
while that in DLAs is closer to the trend for the Sagittarius dwarf. This is 
consistent with the idea that sub-DLAs may have an older stellar population 
that got enriched faster than the stellar population in DLAs. We caution, 
however, that the existing data on [Mn/Fe] in galaxies is limited, and dwarf 
galaxies can have a range of Mn/Fe values depending on their star formation 
history (e.g., Cescutti et al. 2008). 

Consistency with galaxy down-sizing: The star formation history of 
galaxies shows evidence of down-sizing: the more massive galaxies experience 
quenching of star formation activity at earlier times and get enriched faster 
than less massive galaxies (e.g., Cowie et al. 1996; Bundy et al. 2006). If sub- 
DLAs arise in more massive galaxies than DLAs, then the faster metallicity 
evolution of sub-DLAs would be naturally expected from down-sizing. 

Furthermore, the existing sub-DLA metallicity data pose no serious con- 
flict with the observations of massive galaxies being already metal-rich at 
z ~ 2 (e.g., Panter et al. 2008), for the following reasons. 

1. Fig. 1 shows that within the error bars, the mean sub-DLA metal- 
licity at z ~ 2 is consistent with being only 0.6 dex below the solar value. 
Thus, there is no evidence that sub-DLAs evolved much later in metallicity 
compared to massive galaxies. (A suggestion to the contrary in some earlier 
works based on Fe/H vs. z does not account for the fact that the Fe that 
is claimed to be too low could be locked up in dust grains.) Of course, data 
on Zn and S in a larger sample of sub-DLAs would help to address this issue 
further. 

2. Fig. 3 of Panter et al. (2008) shows the metallicity-redshift relation 
for gas in star-forming galaxies. Comparing this figure with our Fig. 1, 
the sub-DLA data at z > 1.1 lie between the green and purple curves in 
the former figure (stellar mass-ranges of 1-3 xlO lo M and < 1 x 10 10 M Q , 
respectively). In other words, existing sub-DLA data are consistent with 
their being moderately massive galaxies. 

3. Erb et al. (2006) showed that there is an offset in the mass-metalicity 
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relations at z ~ and z ~ 2, consistent with our metallicity evolution results, 
within the error bars. 

As noted before, the sub DLAs may not be the most massive galaxies. In 
fact, massive ellipticals may have such hot gas that the absorbers we normally 
know of evaporate, heated by infall; in this case we would not sample such 
galaxies even though they may have super-solar abundances. 

Comparison with metallicity evolution of luminous star-forming 
galaxies: It has been pointed out that star forming galaxies (SFGs) show 
a shallower rise in metalicity with decreasing redshift than do sub-DLAs. 
However, the error on the slope of the SFG relationship is not well enough 
determined to know if that is true or not. The slope of the metallicity 
evolution for SFGs is -0.15 according to Li (2008) , but ps -0.3 at z < 0.4 
according to Lara-Lopez et al. (2009). Thus, it is not clear whether there is 
a significant difference from the slope for sub-DLAs (—0.46 ± 0.18). 

Furthermore, there is no a-priori reason to expect the chemical evolution 
of sub-DLAs to be similar to that of luminous star-forming galaxies. Massive 
(or moderately massive) galaxies are not always forming stars. Sub-DLAs 
could well have already passed through the star formation stage (as expected 
in the down-sizing scenario if they are more massive). Indeed, the few at- 
tempts made so far to detect direct signatures of star formation in absorbers 
at z > 0.7 have often led to non-detections or weak detections of emission 
lines (Ly-a, [O II], H-a, etc.), implying relatively low star formation rates 
(SFRs; e.g., Kulkarni et al. 2000, 2001, 2006, Peroux et al. 2010a, and 
references therein; Rauch et al. 2008; Straka et al. 2010). 

Consistency with Number Densities: The natural, but unfounded, 
association of high H I column density of a single absorbing cloud with the 
entire mass of a galaxy might make the number ratios of sub-DLAs to DLAs 
(a factor of ~ 4) seem to contradict the possibility that sub-DLAs come from 
more massive galaxies. We point out, however, that the number density of 
absorption lines per unit redshift (dN/dz) depends on the product of the 
cross-section of the absorbers, and their number density per unit co-moving 
volume (see, e.g., Sargent et al. 1980). Even though the number density 
of massive galaxies per unit co-moving volume is smaller than that of less 
massive ones, the cross-section for absorption of massive galaxies is higher 
than that of less massive ones. For example, the slope of the galaxy stellar 
mass function is ~ —0.4 in the range 8 < log (M*/M Q ) < 10.3 (e.g., Baldry, 
Glazebrook, & Driver 2008). But the cross-section rises more rapidly with 
mass. Using hierarchical galaxy formation models, Okoshi & Nagashima 
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(2005) show that the average cross-section of an absorber is expected to 
scale as a oc M^ 97 for 6.5 < log (Mhi/M ) < 10.5, which also agrees with 
observations of local H I-rich galaxies (e.g., Rosenberg & Schneider 2003). 
Therefore, the product of galaxy number density and cross-section, and hence 
dN/dz, should increase with the mass. Thus, if sub-DLAs are more massive 
than DLAs, then the fact that sub-DLAs outnumber DLAs could be a result 
of the larger gas cross-section for sub-DLAs and can be consistent with the 
relative fractions of more massive vs. less massive galaxies. Indeed, far more 
sight-lines through the Milky Way would be considered sub-DLAs, than those 
that would qualify as DLAs. 

6. Conclusions and Future Work 

We have examined the metallicity evolution of DLAs and sub-DLAs using 
the largest sample yet of the nearly undepleted elements Zn and S, and 
compared them to properties of galaxies. Our main conclusions are as follows: 

(1) The higher mean metallicities of sub-DLAs compared to DLAs appear 
to be real, and not an artefact of Mg II selection, Fe II selection, or other 
systematic effects. 

(2) The intercept of the metallicity-redshift relation based on Zn or S is 
higher than that based on Fe by 0.5-0.6 dex for both DLAs and sub-DLAs. 

(3) A correlation of metallicity with W 2 796 or W2599 is seen for DLAs, but 
not for sub-DLAs. 

(4) The sub-DLAs have somewhat higher velocity spreads than DLAs, and 
a higher fraction of systems with Av > 150 km s -1 . 

(5) Compared to DLAs, the sub-DLAs show a steeper growth of [Mn/Fe] 
with increasing metallicity that is closer to the trend for Galactic bulge/ 
thick disk stars. 

(6) The absorber data are consistent with galaxy down-sizing. 

(7) The data are also consistent with the relative number densities of low- 
mass and high-mass galaxies. 

(8) It is plausible that sub-DLAs arise in more massive galaxies on average 
than DLAs. 

Why is there a metallicity difference between DLAs and sub-DLAs? We 
suggest the physical reason for this difference is that they differ in their star 
formation histories. Sub-DLAs, being more massive on average than DLAs, 
underwent more rapid star formation and gas consumption, leaving them 
with lower Ahi and higher average metallicity-as also expected from galaxy 
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down-sizing. The difference between sub-DLAs and DLAs may be qualita- 
tively similar to the dichotomy between early-type and late-type galaxies. 

We stress that a commonly hypothesized picture associating absorbers 
of decreasing N m with progressively larger regions around galaxies does not 
agree with the observed metallicity differences between DLAs and sub-DLAs. 
If sub-DLAs were sampling outer parts of a galaxy, one would expect them to 
be less metal-rich than DLAs, in contrast to the observations. In the scenario 
we suggest, sub-DLAs have lower JVhi not because they sample outer parts of 
galaxies, but because they experience more rapid star formation that leaves 
less neutral gas remaining. We also emphasize that there is little evidence 
for the hypothesis associating decreasing N m with larger regions. The exist- 
ing spectroscopically confirmed detections of galaxies producing DLAs/sub- 
DLAs show no systematic increase in impact parameter with decreasing Nm 
(e.g., Peroux et al. 2010a; Meiring et al. 2010). Another way to test this sce- 
nario would be by observing DLAs/sub-DLAs toward gravitationally lensed 
or close binary quasars that sample multiple sight-lines at different impact 
parameters through foreground galaxies. 

Some hints that metal-rich absorbers may arise in massive galaxies are 
seen in the small amount of imaging data available so far. A recent adaptive 
optics imaging study suggests a massive bulge-dominated galaxy as the source 
of a supersolar absorber at z = 0.7 (Chun et al. 2010). Furthermore, based 
on kinematics of H-a emission lines detected in VLT SINFONI integral field 
spectroscopy for two quasar absorbers (one DLA and one sub-DLA), Peroux 
et al. (2010b) estimate the halo mass of the supersolar sub-DLA galaxy to 
be ~ 25 times that of the relatively metal-poor DLA galaxy. Some other 
recent imaging studies have found multiple candidate galaxies near metal- 
rich sub-DLAs (e.g., Hewett & Wild 2007; Straka et al. 2010), suggesting 
that these metal-rich sub-DLAs may be associated with galaxy groups or 
clusters. (We note, however, that in \ow-z galaxies, metallicity shows only a 
modest correlation with environment [e.g., Cooper et al. 2008; Ellison et al. 
2009b]). 

We also note the possibility that the H I column density of a given ab- 
sorber may change with time, as a result of the interaction between stars and 
interstellar gas. For example, star formation and/or gas stripping may cause 
a DLA to lose some H I and turn into a less gas-rich but more metal-rich 
sub-DLA. On the other hand, mergers may cause sub-DLAs poorer in H I to 
evolve into more gas-rich DLAs. Thus, the DLA and sub-DLA populations 
may not always be distinct. 
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Further studies of the number density distribution per unit redshift per 
unit H I column density interval d 2 N / dzdNm at low redshift are essential 
for both DLAs and sub-DLAs, to help understand whether a causal con- 
nection could exist between DLA and sub-DLA populations. Future high 
S/N UV spectroscopy of quasars with low- z sub-DLAs will help toward a 
better understanding of d 2 N/dzdN m by reducing the incompleteness at the 
Iow-A^hi end. Larger abundance samples for nearly undepleted elements like 
Zn and S are essential to better constrain the evolution of sub-DLAs and 
DLAs. Relative abundances such as [Mn/Fe] in a large sample would help to 
better understand the stellar populations in DLAs and sub-DLAs. Finally, 
it is essential to image the galaxies producing the absorbers and measure 
their masses, morphologies, and SFRs to directly understand how DLA and 
sub-DLA populations differ and how they are related to each other. 

Table 1: iy 2 796-metallicity Correlations for DLAs (D) and Sub-DLAs (S). The statistic 
used is the generalized Kendall's r correlation coefficient with survival analysis. 



Type 


Quantity 


Result 


S 


[X/H] a , W 2796 


T = -0.070, P = 0.670 


S 


[Fe/H], W 2796 


r = 0.128, P = 0.541 


D 


[X/H] a , W 2796 


r = 0.819, P < 0.0001 


D 


[Fe/H], W 27m 


r = 0.752, P = 0.010 



X = Zn or S 



Table 2: W2599-metallicity Correlations for DLAs (D) and Sub-DLAs (S). The statistic 
used is the generalized Kendall's r correlation coefficient with survival analysis. 



Type 


Quantity 


Result 




S 


[X/H] a , W 2599 


r = 0.171, P = 


0.302 


D 


[X/R] a , W 2599 


r = 0.744, P = 


0.0005 



X = Zn or S 
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Table 3: Velocity Spread Distributions of DLAs (D) and Sub-DLAs (S) 



Type 


Quantity 


Statistic 


Result 


S 


Av 


(weak lines) 


Mean, Median 


125.5, 95.3 km s" 1 


D 


Av 


(weak lines) 


Mean, Median 


103.0, 80.0 km s- 1 


S 


Av 


(Mg II 2796) 


Mean, Median 


225.0, 197.0 km s" 1 


D 


Av 


(Mg II 2796) 


Mean, Median 


191.8, 177.0 km s" 1 


S 


Av 


(weak lines) 


Fraction with 


30.4 ±8.1% 








Av > 150 




S 


Av 


(weak lines) 


Fraction with 


17.4 ±6.1% 








Av < 50 




D 


Av 


(weak lines) 


Fraction with 


11.9 ±4.2% 








Av > 150 




D 


Av 


(weak lines) 


Fraction with 


29.9 ± 6.7% 








Av < 50 
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Table 4: Velocity-Metallicity Correlations for DLAs (D) and Sub-DLAs (S). The statistic 
used is the generalized Kendall's r correlation coefficient with survival analysis. 



Type 


Quantity 


Result 


D 


[X/HT, Av (weak) 

L / J J \ / 


r = 0.914, 






P < 0.0001 


S 


[X/H]°, Av (weak) 


r = 0.106, 






P = 0.526 


D 


[Fe/H], Av (weak) 


r = 0.960, 






P < 0.0001 


S 


[Fe/H], Av (weak) 


r = 0.498, 






P = 0.0277 


D 


[X/H], Av (Mg II) 


r = 0.765, 






P = 0.017 


S 


[X/H], Av (Mg II) 


r = -0.097, 






P = 0.638 



X = Zn or S 
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Figure 1: A(HI)-weighted mean metallicity vs. look-back time relation for 154 DLAs and 
58 sub-DLAs with Zn or S measurements. Filled circles show 8 bins with 19 or 20 DLAs 
each. Unfilled circles refer to the lowest time bin for DLAs split into 2 bins with 10 DLAs 
each. Squares denote 4 bins with 14 or 15 sub-DLAs each. Horizontal bars denote ranges 
in look-back times covered by each bin. Vertical errorbars denote 1 a uncertainties. The 
triangle shows the formal lower limit to the average metallicity for a composite spectrum 
from 698 absorbers with average log Ahi ~ 20 (sample 24) from York et al. (2006). The 
bold solid and dashed curves show the best fits obtained from linear regression of the 
metallicity vs. redshift data for sub-DLAs and DLAs, respectively. The light dot-dashed 
and dot-double-dashed curves show, respectively, the mean metallicity in the models of 
Pei et al. (1999) and Somerville et al. (2001). Sub-DLAs appear to be more metal-rich 
and faster-evolving than DLAs, especially at lower redshifts. 



23 



.5 -2 -1.5 -1 -0.5 0.5 

[Zn/H] 



Figure 2: [S/Zn] vs. [Zn/H] for absorbers with detections of S and Zn. Dashed horizontal 
lines denote the solar level and ±0.2 dex levels. S/Zn tracks the solar ratio within ±0.2 
dex for most of the absorbers. 
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Figure 3: Rest-frame equivalent width of Mg II A2796 plotted vs. metallicfty [X/H] from 
Zn or S, for (a) sub-DLAs (top) and (b) DLAs (bottom). No correlation is seen for sub- 
DLAs, whereas a correlation is seen for DLAs. Thus sub-DLA data do not have an Mg 
Il-selection bias toward high-metallicities. 
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Figure 4: Comparison of data and simulations for metallicity vs. rest-frame equivalent 
width of Fe II A2599, for (a) sub-DLAs {top), and (b) DLAs (bottom). The squares, 
triangles, and circles show W2599 plotted vs. metallicity [X/H] for the absorbers in our 
sample where W2599 data are available. The diamonds show median W2599 values from 
30,000 simulations, plotted as a function of metallicity, in bins of 0.25 dex in [X/H]. The 
dashed horizontal error bars on the diamonds show the ±1<7 ranges in the W2599 values 
generated in the simulations. The DLA data show a correlation between the metallicity 
and W2599, but the sub-DLA data do not show a correlation. Thus, sub-DLA data do not 
have an Fe II selection bias toward high metallicities. 
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Figure 5: Velocity spread Av, determined from weak absorption lines, plotted vs. metal- 
licity [X/H] from Zn or S, for (a) sub-DLAs and (b) DLAs. The correlation is stronger for 
DLAs than for sub-DLAs. 
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